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Abstract 

66" 
°V 

> \ , We use equations for Faddeev amplitudes to solve the shell-model problem for 

four nucleons in the model space that includes up to lAhQ harmonic-oscillator 
| excitations above the unperturbed ground state. Two- and three-body effec- 

Q ■ tive interactions derived from the Reid93 and Argonne V8' nucleon-nucleon 

potentials are used in the calculations. Binding energies, excitations energies, 
CN , point-nucleon radii and electromagnetic and strangeness charge form factors 

for He are studied. The structure of the Faddeev-like equations is discussed 
J> | and a formula for matrix elements of the permutation operators in a harmonic- 

oscillator basis is given. The dependence on harmonic-oscillator excitations 



(N 



■ allowed in the model space and on the harmonic-oscillator frequency is investi- 

gated. It is demonstrated that the use of the three-body effective interactions 
00 '. improves the convergence of the results. 



O 

I. INTRODUCTION 



X 



Many different methods have been used to solve the few-body problem in the past. One 



of the most viable approaches appears to be the Faddeev method pL It has been suc- 



cessfully applied to solve the three-nucleon bound-state problem for various nucleon-nucleon 
potentials . For solution of the four-nucleon problem one can employ Yakubovsky's 

generalization of the Faddeev formalism J7j as done, e.g., in Ref. ||. Alternatively, other 
methods have also been succesfully used in the past, such as, the Green's function Monte 
Carlo method or the correlated hyperspherical harmonics expansion method [|10|| . 

On the other hand, when studying the properties of more complex nuclei one typically 
resorts to the shell model. In that approach, the harmonic-oscillator basis is used in a 
truncated model space. Instead of the free nucleon-nucleon potential, one utilizes effective 
interactions appropriate for the truncated model space. Examples of such calculations are 
the large-basis no-core shell-model calculations that have recently been performed |Tl|Jl~^1 . In 
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these calculations all nucleons are active, which simplifies the effective interaction as no hole 
states are present. The effective interaction is determined for a system of two nucleons in 
a harmonic-oscillator well interacting by the nucleon-nucleon potential and is subsequently 
used in the many-particle calculations. 

In a recent paper we combined the shell-model approach to the three-nucleon problem 



with the Faddeev formalism [ 13| . That allowed us to extend the shell-model calculations to 
a model space of excitations of 32HQ above the unperturbed ground state and to study the 
convergence with respect to the size of the model space. In the present paper we generalize 
these earlier calculations to the four-nucleon problem. We introduce equations for Faddeev 
amplitudes that are fully antisymmetrized for three nucleons. As the center-of-mass term 
is removed, we are able to work in a model space up to an excitation of 14hQ above the 
unperturbed ground state. For comparison the largest shell-model calculations so far for 
4 He are those performed by R. Ceuleneer et al., in which a 10hQ model space was utilized 
rj|. The main motivation for the present work is to test the shell-model approach and the 



effective interactions derived from realistic nucleon-nucleon (NN) potentials that are used in 
the conventional shell-model applications for more complex systems. As the equations that 
we employ can be conveniently used with three-body interactions or three-body effective 
interactions, we investigate, in addition to two-body effective interactions, also three-body 
effective interactions in the present formalism. Such effective interactions are not typically 
used in the traditional applications. We show that the inclusion of three-body effective 
interactions improves the overall convergence of the results. At the same time our work serves 
as an alternative method to solving the four-nucleon problem. We can study the convergence 
properties of the results with the increasing size of the model space. If convergence is 
achieved, our results will approach the exact solutions obtained by other methods. In our 
formalism we seek simultaneously solutions for both the ground-state and the excited states. 
In the past, the variational Monte Carlo method was used to investigate the excited states 
of 4 He using realistic NN potentials WE\. In most four-nucleon calculations with realistic NN 



potentials, however, only the ground-state properties were evaluated [jq, |10|1 . On the other 
hand, earlier studies that investigated the excited-state properties usually did not employ 
realistic NN potentials |H|,IB[. Recently, the four-nucleon resonant and scattering states 



were investigated using realistic NN potentials in the framework of the resonating group 
method [H], the correlated- hyperspherical-harmonics method as well as in the solution 



of the Faddeev-Yakubovsky equations in configuration space . 

The present calculation is simplified by using a compact formula for the matrix elements 
of the permutation operators in the harmonic-oscillator (HO) basis. Also, because of the way 
we do the model-space truncation, we keep equivalence of the Faddeev-like and Schrodinger 
equations throughout the calculation. In addition to calculation of ground-state and ex- 
cited state energies and point-nucleon rms radii, we also evaluate electromagnetic (EM) and 
strangeness form factors in the impulse approximation. 

In section |I| we first discuss the Faddeev equations for the shell-model problem of three 
nucleons. Then a generalization to the four-nucleon system is introduced. In section [TTI| we 
present the energy, radii and form factor results for 4 He. Conclusions are given in section 

m 
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II. SHELL MODEL AND FADDEEV-LIKE FORMALISM 

In shell-model studies the one- plus two-body Hamiltonian for the A-nucleon system, 

i.e., 

# = E 4 + - 3) , (1) 

i=i zm i<j 

where m is the nucleon mass and V^(fi — fj), the NN interaction, is usually modified by 
adding the center-of-mass HO potential ^AmO, 2 R 2 , R = -jJ2t=i^i- This potential does 
not influence intrinsic properties of the many-body system. It provides, however, a mean 
field felt by each nucleon and allows us to work with a convenient HO basis. The modified 
Hamiltonian, depending on the HO frequency Q, can be cast into the form 

A 

1=1 

The one-body term of the Hamiltonian (0) can be re- written as a sum of the center-of-mass 

term, H^ m = + \AmVt 2 R 2 , where P cm = J2f=iPi, an d a term depending only on the 
relative coordinates. In the present application we use a basis, which explicitly separates 
center-of-mass and relative-coordinate wave functions. Therefore, the contribution of the 
center-of-mass term is trivial and will be omitted from now on. 

The shell-model calculations are performed in a finite model space. Therefore, the in- 
teraction term in Eq. (|2|) must be replaced by an effective interaction. In general, for an 
A-nucleon system, an A-body effective interaction is needed. In practice, the effective in- 
teraction is usually approximated by a two-body effective interaction. In the present study 
we will also employ a three-body effective interaction. As approximations are involved in 
the effective interaction treatment, large model spaces are desirable. In that case, the cal- 
culation should be less affected by any imprecision of the effective interaction. The same is 
true for the evaluation of any observable characterized by an operator. In the model space, 
renormalized effective operators are required. The larger the model space, the less renormal- 
ization is needed. We may take advantage of the present approach to perform shell-model 
calculations in significantly larger model spaces than are possible in the conventional shell- 
model approach. At the same time we can investigate convergence properties of effective 
interactions. 
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A. Three-nucleon system 



In this subsection we repeat the steps discussed in Ref. |TB 
three-nucleon shell-model problem in the Faddeev formalism 
i.e., A = 3, the following transformation of the coordinates 



that are needed to solve the 
For a three-nucleon system, 



and, similarly, of the momenta, can be introduced that brings the relative-coordinate part 
of the one-body HO Hamiltonian into the form 

H = ^ + -m^V + ^- + -mQ 2 f . (4) 
2m 2 2m 2 y V ; 

Eigenstates of this Hamiltonian, 

\nlsjt,MCJ,JT) , (5) 

are then used as the basis for the three-nucleon calculation. Here n, I and A/", C are the HO 
quantum numbers corresponding to the harmonic oscillators associated with the coordinates 
and momenta r,p and y,q, respectively. The quantum numbers s,t,j describe the spin, 
isospin and angular momentum of the relative-coordinate partial channel of particles 1 and 
2, while J is the angular momentum of the third particle relative to the center of mass 
of particles 1 and 2. The J and T are the total angular momentum and the total isospin, 
respectively. 

The Faddeev equation for the bound system can be written in the form 

H\<j>) = E\cf>) , (6) 

with 

H = H + V(r)T . (7) 
Here, V(r) = V N (y/2r) - \mQ 2 r 2 is the potential and T is given by 

r = i + r(-> + r (+) , (8) 

with 7~( + ) and 7~' ' the cyclic and the anticyclic permutation operators, respectively. Pre- 
viously jnj, we derived a simple formula for the matrix elements of + in the basis 
(|), namely 

{nxhsxjxbxMxdJi, JT\T { - ] + T {+) \n 2 l 2 s 2 j 2 t 2) M 2 C 2 J 2) JT) = -6 Nl , Na 

\ h sx jx 

xY,L 2 S 2 ^ 2 JxJ 2 sxS 2 ixi 2 {-l) L \Cx \ Jx 

ls {LSJ){LSJ)^ 2 
_ 1)sl+S2+tl+i2 _ £l _^ ^^^^^i^^^^a + (nxZiM^i^lA^rsna^a^a] , (9) 

where N 4 = 2ui + k + 2J\f { + £ h i = 1,2, j = y/2j + 1 and {MiCxn l lxL\n 2 l 2 N 2 £ 2 L) ? , is the 
general HO bracket for two particles with mass ratio 3 as defined, e.g., in Ref. ||20|| . The 
expression (|9]) can be derived by examining the action of and T^~' on the basis states 
(H). A similar derivation for a different basis is described, e.g., in Refs. JH],^]. Let us note 




that it follows from the antisymmetry of the two-nucleon states and from the symmetry 
properties of the HO brackets that the contributions of and in (§) are identical. 

The eigensystem of the operator T (W) consists of two subspaces. The first subspace 
has eigenstates with the eigenvalue 3, which form totally antisymmetric physical states. 
The second subspace has eigenstates with the eigenvalue 0, which form a not completely 
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antisymmetric, unphysical subspace of states. We found these properties of T by direct 
calculation using the relation (||). It is, however, a general result. The same structure of 



eigenstates was also obtained in Ref. [23|] using a different basis. The eigenvalue structure 
follows from the fact that |T has the properies of a projection operator. It is possible to 
hermitize the Hamiltonian (|7|) on the physical subspace, where it is quasi-Hermitian. The 
Hermitized Hamiltonian takes the form 

H = H + T 1/2 V(f)f 1/2 , (10) 

where T operates on the physical subspace only. 

The operator T (H) is diagonal in N = 2n + I + 2J\f + C. Note that any basis truncation 
other than one of the type N < iV max violates, in general, the Pauli principle and mixes 
physical and unphysical states. Here, iV max characterizes the maximum of total allowed HO 
quanta in the model space and is an input parameter of the calculation. The truncation 
into totally allowed oscillator quanta iV < iV max , however, preserves the equivalence of the 
Hamiltonians (|7]) and ( [TD|) on the physical subspace. 



B. Four-nucleon system 

By relying on the results obtained for the three-nucleon system, as described in the 
previous subsection, we can extend the formalism to the four-nucleon system. We use the 
Hamiltonian (H) with A = 4. By introducing the coordinate (and momentum) transforma- 
tions, 



yhin - r 2 ) , (11a) 



2 

V= ^HW^+^-rz], (lib) 
z= ^Wi +r 2 + f 3 ) -r 4 ] , (11c) 

we obtain the one-body part of the Hamiltonian (Q) in the form 

H = ^ + -mQ 2 ^ + ^- + -mQ 2 f + — + -mQ 2 ^ , (12) 
2m 2 2m 2 y 2m 2 ' K J 

with the center-of-mass term omitted. 

A possible generalization of the Faddeev equation (0) for four identical particles can be 
written in the form 

##(123)4) = £#(123)4) , (13) 

with 

##(123)4) = 1^(123)4) + i(Vl2 + Via + V 23 ) (1^(123)4) + |'0(432)l) + 1^(134)2) + 1^(142)3)) , 

(14) 

and 
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(|-0(123)4> + |^(432)l) + 1^(134)2) + 1^(142)3)) = (1 - Tia - T 24 ~ T34) |V(123)4> = %\^{12Z)i) • 

(15) 

Here, 1^(123)4) is a four-fermion Faddeev amplitude completely antisymmetrized for particles 
1,2, and 3. There are three other equations that can be obtained from Eq. ([13]) by permut- 
ing particle 4 with particles 1, 2, and 3. Their sum then leads to the Schrodinger equation. 
We note that the present equations are different from the traditional Faddeev-Yakubovsky 
equations [0], which combine Faddeev amplitudes depending on two sets of relative coor- 
dinates. We are working with a complete orthonormal basis. It is, therefore, sufficient 
and convenient to use a single set of coordinates defined by the relations (0). Unlike the 
Faddeev amplitudes used typically in the Faddeev-Yakubovsky equations, the amplitudes 
appearing in Eq. ( |i~3|) are antisymmetrized with respect to the first three particles. Those 
amplitudes are obtained, as described below, in a straightforward manner with the help 
of our three-nucleon HO formalism introduced earlier. The present equations allow us to 
employ easily real three-body interactions or three-body effective interactions. The latter 
property makes them particularly useful for the present extension of shell-model calculations 
for four nucleons. At the same time, the use of Faddeev amplitudes antisymmetrized for 
particles 1, 2 and 3 allows us to reduce the dimmension of the basis significantly. 
We start the four-nucleon calculation using the basis 

imiJtfunJiJt, JT) . (16) 

with the three-fermion part given by the antisymmetrized eigenstates of T (§) corresponding 
to eigenvalue 3, e.g., 

\N l iJ l T l )=Y,cXmcjMlsjt,M^JiTi) , (17) 

where N± = 2n + / + 27V + C and i counts the eigenstates of T with the eigenvalue 3 for 
given Ni and Ji,T±. Further, n z ,l z are the HO quantum numbers corresponding to the 
harmonic oscillator associated with the coordinate z and the momentum o and J± is the 
angular momentum of the fourth particle relative to the center of mass of particles 1, 2 and 
3. 

As in the case of the three-particle transposition operators (0), a compact formula can 
be derived for the matrix elements of the four-particle transposition operators in the basis 
(0),e.g., 



(^lL^L^lL^lL, n z\}zhJ^-, JT\T 14: + 7^4 + % 4 \ -/V 1R 2 R J 1R T 1R , n 2R / 2Rl 74 R , JT) 

— A N 1L i L J 1L T 1L N m i R J 1R T 1R f 2 jl 

- °A f L,A r R Z^ C nLiLSLiLtLA^£L^3L C ^R5RiRtRAr R £ Rt 7 3H ^lL^l 

XjLjRi3L l 73Ri4L l 74RJlLilRT 1L f 1R (-l) Tl ^ TlR+5lL+5lR 
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ir R+ilL (-l)^(n / rn, L Z, L L> zR UMii:R^)8(riLiLM J >C LJ L 1L |n / Z / n R i R L 1L ) 
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-^R^L,SR^L,^A/^^^ > (18) 

where /V"x = 2rix + ^x + 2A/x + £x + 2n 2 x + 4x,X = L or R, and, e.g., the expression 
{n z \Jz\J^^\,L'\M^.C^n z -^l z ^,L')^ denotes a general HO bracket for two particles with mass 
ratio 8, as defined in Ref. [p0| |. Similarly, as in Eq. (|9]) the brackets for two particles with mass 



ratio 3 also appear in the relation (jl8]) . In the derivation of the expression we relied on 
the antisymmetry of the basis states with respect to particles 1,2 and 3. The calculation was 
facilitated by application of the operators — 7^ and — 7^ 3 . The relation ( |i~8]) appears to be 
non-symmetric. However, its numerical evaluation leads to a symmetric matrix. It may also 
appear that the angular momentum sums in (|18"D can be summed up. In fact, it is possible 
to simplify the expression by introducing a 15j-coefficient of the fifth kind as defined, e.g., in 
Ref. |[24|| , but, as such coefficients are seldomly used, we prefer to keep the summations in the 



explicit form. On the other hand, a significant simplification of the expression ([18]) can be 
obtained, when the symmetry relations of different terms are exploited. First, it follows from 
the properties of the HO brackets and from the antisymmetry of the two-nucleon states that 
the contributions of 7i 4 and are identical. Second, using the fact that the states ([16]) are 
antisymmetrized for the particles 1, 2 and 3 it follows that all three permutation operators 
appearing in ( |T8"D give identical contributions to the expression (JT8|) . The computation of 
7J34 is the simplest. In that case a partial summation of the angular momentum coefficients 
can be performed, yielding a compact expression 



(-/ViL^L JixTlL, ^zL^lJ^l, JT\T^\Nmi^ JirTir, UzrIzrJaRi JT) 
- °A r L,A r R 2^ C nlsjWi J C^J3i J C nlsjW R C B .J3B. 
xi 3 Li3Ri4LJ4RJlLJlRT 1L r iR (-l) TlL+TlR+ ^ L+ ^ 



/ £ L i zR K \ j C R l z 



j J3L JlL 

xK 2 { J m K J n . 
xL /2 (-l)^ + ^ +i \n zL ^ L £ L Lq^ R £ R n, R U^)8- (19) 



1 
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t 




) 




1 
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T 


Til 






K 




\ /zL 
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\ ^zR 


£l 


V 



Thus, by multiplying the expression (|T^) by three we obtain the same matrix element as from 
(|18|). We note that a generalization of the evaluation of the permutation operator matrix 



element (19) to a more complex system, than the presently studied A = 4, is straightforward. 
Its simplicity suggests that the present formalism can be extended to systems with A > 4. 
Similarly, as for the operator T (H), eigenstates of the operator T4 defined by the relation 



(|i5|) can be subdivided into two subspaces. A physical subspace is spanned by totally 
antisymmetric states, in this case corresponding to the eigenvalue 4, and a spurious subspace 
is spanned by eigenvectors corresponding to the eigenvalue 0. It is possible to symmetrize 
the Hamiltonian H appearing in Eq. ( |r3"P on the physical subspace. The symmetrized 
Hamiltonian then takes the form 

H = H + % 1,2 \{V 12 + Vi 3 + V 23 )% 1/2 , (20) 
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where 7^ operates only on the physical subspace. In our calculations, described later, we 
diagonalize the symmetrized Hamiltonian ( p0|) in the physical basis formed by the eigenstates 
ofT 4 . 

The operator 7^ ( |i"5|) is diagonal in N = 2n + I + 2j\f + C + 2n z + l z . A basis truncation 
defined by a restriction on the totally allowed oscillator quanta N < iV max preserves the 
equivalence of the Hamiltonians (^) and (pTOf) on the physical subspace. 



C. Effective interactions 

From solving two-nucleon systems in a HO well, interacting by soft-core potentials, one 
learns that excitations up to about 300HQ (iV max = 300) are required to get almost exact 
solutions. We anticipate, therefore, that at least the same number of excitations should be 
allowed to solve the many-nucleon system. The Faddeev formulation has the obvious advan- 
tage compared with the traditional shell-model approach that the center-of-mass coordinate 
is explicitly removed. Even then, it is not feasible to solve the eigenvalue problem either 
for ( |TUD or for fl2T)D in such a large space. On the other hand, shell-model calculations are 
always performed by employing effective interactions tailored to a specific model space. In 
practice, these effective interactions can never be calculated exactly, because, in general, for 
an A-body effective interaction is required for an A-nucleon system. We may, however, take 
advantage of the present approach to perform shell-model calculations in significantly larger 
model spaces than are possible in conventional shell-model approach. At the same time we 
can investigate convergence properties of effective interactions. If convergence is achieved, 
we should obtain the exact solution, since by construction the effective interactions that we 
employ satisfy the condition V e g — > V for iV max — ► oo. 

Usually, the effective interaction is approximated by a two-body effective interaction de- 
termined from a two-nucleon system. In the present calculations we replace matrix elements 
of the potential V(r) by matrix elements of an effective two-body interaction, derived in 
a straightforward manner for each relative-coordinate partial channel. The relevant two- 
nucleon Hamiltonian is then 

H 2 = H 02 + V = + \mtff k + V N (V2r) - ^f 2 , (21) 

which can be solved as a differential equation or, alternatively, can be diagonalized in a 
sufficiently large harmonic oscillator basis. For a four-nucleon system we set A = 4 in Eq. 
(pl|), which implies that we are dealing with a bound-state problem. 

To construct the two-body effective interaction we employ the Lee-Suzuki |25j similarity 



transformation method, which gives the effective interaction in the form PV c rP = PVP + 
PVQujP, with uj the transformation operator satisfying u = QujP, and P and Q = 1— P, the 
projectors on the model and the complementary spaces, respectively. Our calculations start 
with exact solutions of the Hamiltonian (|2l|) and, consequently, we construct the operator uj 
and, then, the effective interaction directly from these solutions. Let us denote the relative- 
coordinate two-nucleon HO states, which form the model space, as \ap), and those which 
belong to the Q-space, as |«q). Then the Q-space components of the eigenvector \k) of the 
Hamiltonian ([H]) can be expressed as a combination of the P-space components with the 
help of the operator uj 
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(a Q \k) = ^(a Q \uj\ap)(ap\k) . (22) 

Ctp 

If the dimension of the model space is dp, we may choose a set K, of dp eigenevectors, 
for which the relation (p2|) will be satisfied. Under the condition that the dp x dp matrix 
(ap\k) for \k) G K, is invertible, the operator u can be determined from (0). In the present 
application we select the lowest states obtained in each channel. Their number is given by 
the number of basis states satisfying In + 1 < iV max . Once the operator uj is determined, the 
effective hamiltonian can be constructed as follows 



(jp\H 2eS \a P ) = 



fce/c 



(j P \k)E k (k\ap) + ^2(~fp\k)E k (k\a Q )(aQ\Lj\ap) 



(23) 



It should be noted that P\k) = J2 ap \ a p)( a p\k) for \k) G /C is a right eigenvector of ( |23|) 
with the eigenvalue E k . 

This Hamiltonian, when diagonalized in a model-space basis, reproduces exactly the set 
/C of dp eigenvalues E k . Note that the effective Hamiltonian is, in general, quasi-Hermitian. 
It can be hermitized by a similarity transformation determined from the metric operator 
P(l + uj^uj)P. The Hermitian Hamiltonian is then given by 



H- 



2cS 



P(l + Ju)P 1/2 H 2cS P(l +Ju)P ^\ (24) 



-1/2 



Finally, the two-body effective interaction used in the present calculations is determined 
from the two-nucleon effective Hamiltonian fl24"|) as V 2e fi = H 2cS — H 02 . We note that the 
interaction V\ 2 + + V 23 in Eq. ( p0|) is then replaced by T 1 / 2 V 2c qT 1 / 2 , which is evaluated 
in a straightforward way in the basis (|T6|). 

As pointed out before, the structure of the Hamiltonian (p0|) allows us to employ easily 
three-body effective interactions in addition to the above discussed two-body effective inter- 
actions. We can replace V\ 2 + Vi 3 + V 23 in Eq. (|20|) by V3 C fj that can be derived from the 
three-nucleon solutions in a similar manner as the two-body effective interaction is derived 
from the two-nucleon solutions. To find we solve the three-nucleon system described 
by the Hamiltonian (0) with V(r) = V N (y/2f) - ^mtfr 2 . As A = 4 we are dealing with a 
bound three-nucleon problem. It can be solved in a three-nucleon model space characterized 
by iV3 max « 30 [|13|]. First, we compute the two-body effective interaction appropriate for the 
model space defined by iV 3max , as discussed earlier in this subsection. Then the three-nucleon 
system is solved in the same space. Afterwards we construct the three-body effective inter- 
action for a model space defined by iV max < iV 3mas . In the present paper we use model spaces 
up to iV max = 14. The effective interaction is constructed exactly, as described above, using 
Eqs. (H-H-EID with H 2eS replaced by H 3eS . The energies E k and the states \k) correspond 



to the three-nucleon system eigenstates, however, and the states \ap) and |q;q) are three- 
nucleon basis states ([[7]) with the model-space condition N\ = 2n + I + 2Af + C < iV max 
and the Q-space condition iV max < N\ < N 3max . The three-body effective interaction is 
computed for different three-nucleon channels characterized by J\ , T\ and parity and is ob- 
tained from the hermitized effective Hamiltonian as V^s = H^es — Ho, where Hq is given 
by Eq. (^). The interaction V^g then replaces V\ 2 + V13 + V23 in Eq. ( pCf) . We note that 
by construction in the limit A^ max — > A^ 3max the three-body effective interaction approaches 
the two-body effective interaction V 3e g — > T l l 2 V 2e tfT 1 / 2 and with A^" max — > 00 the effective 
interaction approaches the bare interaction V 2c q — > V. 
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III. APPLICATION TO 4 HE 



In the present paper we use the Reid93 NN potential |^7| and the Argonne V8' NN 
potential, introduced in Ref. 0. We work in the isospin formalism; the charge invariant 
potential Vn = \V PP + \V nn + \V np is used for each T = 1 wave in the calculations with the 
Reid93 potential. The Coulomb potential is added to V pp in this case. On the other hand, 
the calculations with the Argonne V8' potential, which is isopin invariant, do not include 
the Coulomb potential. 

Our calculation progresses in several steps. The model space is characterized by the 
condition N < iV max , N = 2n+l+2Af+£+2n z +l z . First, the three-nucleon antisymmetrized 
basis is constructed by diagonalizing T (H) in the basis (§) for all Ni = 2n+l+2Af+£ < iV max 
and all J\ , Ti . Then the four-nucleon antisymmetrized basis is calculated by diagonalizing 
T4 (|I5D in the basis (|T6|) for N = Ni + 2n z + l z < iV max with N even for positive parity states 
and N odd for negative parity states. We present results for J = and T = only, but for 
both parities. We note that the four-nucleon basis computation is independent of Q and is 
done only once. The next step is the effective interaction calculation. The two-body effective 
interaction is derived from the Eqs. (p2j)-(|24|). The condition for the relative-coordinate two- 
body effective-interaction model space is then 2n + l < iV max . When solving the two-nucleon 
relative-coordinate Hamiltonian ([H]) in the full space, we truncate the HO basis by keeping 
the states typically up to n = 152. The two-body effective interaction is constructed for all 
partial- wave channels up to j = 6. The resulting effective interaction is finally used as input 
for the four-nucleon calculation, where the Hamiltonian ( p0|) is diagonalized. Instead of a 
two-body effective interaction, we may use a three-body effective interaction, as discussed 
in the previous section. The three-body effective is computed only for the most important 
three-nucleon channels J{F\. In particular, we evaluated the three-body effective interaction 
for J\ = 1/2, 3/2, T\ = 1/2 and for both positive and negative parity. For the channels with 
higher Ji the two-body effective interaction corresponding to iV max is used instead. For 
the parameter iV3 max characterizing the three-nucleon full space, we used A^ 3max = 28 for 
Ji = 1/2 and A 3max = 24 for Ji = 3/2. We also performed calculations with the inclusion 
of the three-body effective interaction for J\ = 5/2 and found it to have little effect. 

Let us remark that the present method for solving the four-nucleon shell-model problem 
is fully equivalent to the standard shell-model approach. In particular, it is straightfor- 
ward to transform the relative-coordinate two-body effective interaction used in the present 
calculations to the two-particle basis used for the shell-model input by the standard trans- 
formation [p8]| . We used the transformed interactions for the model spaces up to an 8hQ 



space to test our results. The shell-model diagonalization was then performed by employing 



the Many-Fermion-Dynamics Shell-Model Code ||29|| , which can be utilized for calculations 



with model spaces comprising up to 9 major HO shells, i.e., iV max = 8 for 4 He. We obtained 
the same results from both the Faddev-like calculation and the standard shell-model cal- 
culation. The Faddeev-like calculation has, obviously, much smaller dimension and can be 
extended to larger model spaces. We also note that we applied the discussed formalism to 
four-electron system in a related study recently pOfl . Our results compared well with those 
obtained by the Stochastic Variational Method j3l| . 
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A. Energies and point-nucleon rms radii 



Our results for the ground-state and excited-state energies and point-nucleon rms radii 
are presented in Figs. [L|-|| where the dependencies on the model-space size and the HO 
energy are shown. A summary of the largest model-space (A max = 14 for the positive- 
parity states and iV max = 13 for the negative-parity states) results is given in Table 0. Let 
us mention an unusual feature of the present calculations, namely, the convergence from 
below for the ground-state energy. It is caused by the asymmetric treatment of the HO 
terms that are added and subtracted to the Hamiltonian in the process of evaluating the 
effective interaction. Our effective interaction is computed for a two- or three-nucleon system 
bound in an HO potential. Therefore, artificial binding from this potential is included in the 
effective interaction and the four-body effects coming from the entire four-nucleon calculation 
may not completely compensate for this spurious binding in a particular model space. We 
note that this type of over-binding in the no-core shell-model calculations was noticed in 
previous studies [^,0,|1^]. This effect decreases as the model-space size increases, as is 
demonstrated in our earlier three-nucleon shell-model calculations [13|. 

In Fig. [I] we present the calculated dependence of the ground-state energy and the 
first-excited + state energy on the model-space size, characterized by A max . The two- 
body effective interaction employed was derived from the Reid93 NN potential. Results for 
hfl = 14, 17, 19, 22 MeV are shown. The corresponding dependence of the point-nucleon rms 
radius is presented in Fig. |^. A slow convergence with the increasing model-space size can 
be observed for the energies with a significantly faster rate for the ground state compared to 
the first excited + state. Also, much stronger dependence of the excited state on the HO 
energy HQ is apparent. The results of the point-nucleon rms-radius calculation demonstrate 
even more the differences between the ground state and the first excited + state. While the 
ground-state radius has almost converged and shows little HQ dependence, the first excited 
+ state displays a strong dependence of its energy on HQ and a steady increase of its 
radius with increasing model-space size. 

Let us remark that in our approach we obtain the ground state as well as the excited 
states by diagonalizing the Hamiltonian. This implies that the excited states are expanded 
in the same harmonic-oscillator basis used for the ground state. While such an approach 
has technical advantages, it might not be physically sound. Cautious interpretation of the 
excited-state results is, therefore, necessary. The significantly different convergence rate of 
the ground state and of the first excited + state manifests the different nature of the two 
states. Let us note that if the model-space size was increased up to the point at which total 
convergence of the excited state was achieved, our procedure would yield isolated three- 
and one-body clusters with an infinite rms radius and a total energy of the three-nucleon 
system. It is possible, though, that we could observe a meta-stability prior to the onset of 
the cluster separation, as the resonance is sharp and low-lying. The present model-space 
sizes, however, are not yet sufficient to arrive at that point. That we have not reached this 
point can be seen from the lack of convergence and, in particular, from the rather small rms 
radius, which shows a significant increase with iV max and a strong dependence on Q. 

The importance of the three-body effective interaction can be judged from the results 
shown in Fig. |^. The ground-state and excited + state energies obtained in a calculation 
that employs the three-body effective interaction is compared to a calculation performed 
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by using only the two-body effective interaction. Results for two different values of the 
HO energy, HQ = 17 and 19 MeV, are presented. The dashed lines connect the two-body 
effective interaction calculation results that are identical to those in Fig. |1] that correspond 
to HO energies of hQ = 17 and 19 MeV. The full lines connect the results obtained in 
calculations with the three-body effective interaction. It is apparent that the three-body 
effective interaction improves the converegence considerably. It is especially true for the 
ground state. The difference between the iV max = 6 and iV max = 14 energies is significantly 
smaller in the calculation that employs the three-body effective interaction. It can also 
be seen that the two-body effective interaction results approach the three-body effective 
interaction results in the largest spaces used in our calculations. In addition, the dependence 
on the HO energy decreases in the three-body effective interaction calculation compared to 
the two-body effective interaction calculation. This holds for both the ground state and the 
first excited + state. However, the inclusion of the three-body effective interaction clearly 
has a larger overall impact on the ground-state results. 

The influence of the three-body effective interaction on the point-nucleon rms radius is 
depicted in Fig. f|. Again we observe a better stability of the radii computed using the three- 
body effective interaction. In particular, the ground-state point-nucleon rms radius shows 
convergence in both the model-space-size dependence and the HO-frequency dependence. 
On the other hand, the three-body effective interaction does not improve the convergence 
of the excited state in any significant way in the model spaces that we employed. 

In Fig. H we show the calculated energies of the first 0~0 state obtained using two- 
body effective interactions in model spaces up to iV max = 13. For a comparison, the results 
for the ground state and the first excited + state from Fig. [3] are also presented. It is 
interesting to note that the 0~0 state shows a better convergence and stability with respect 
to the iV max change as well as a weaker dependence on hfl than the first excited + state. 
This observation is confirmed also in the point-nucleon rms radius calculation as can be 
seen in Fig. || In the experiment, the 0~0 excitation energy, 21.01 MeV, is higher than 
the excitation energy of the first + state, 20.21 MeV. Though in our calculations their 
positions are reversed, it is visible from Fig. || that the extrapolation to larger iV max leads 
to correct ordering of the two states. A possible interpretation of this observation is that 
the excited + state is associated with a radial excitation and, thus, it is more sensitive to 
the HO basis used in our calculations. 

The energy and radius results, obtained using the Argonne V8' NN potential, are pre- 
sented in Figs. and respectively. The three-body effective interaction was used in 
calculating these results, for three different HO energies, Ml = 16, 19 and 22 MeV. The dot- 
ted line represents the value -25.92 MeV obtained for the ground state, using the GFMC pi 



Similarly, as in the calculations with the Reid93 NN potential, we get the best convergence 
for the ground state for the highest value of HQ, while for the excited state the best results 
are obtained for the lowest HQ. The same discussion, given earlier, for the excited-state 
convergence using the Reid93 NN potential, is also valid for the calculations using the V8' 
NN potential. The energy convergence is very slow and there is no sign of convergence of the 
point-nucleon rms radius of the excited + state. A significant dependence on hQ prevails 
for all the model-spaces studied. On the other hand, the ground-state energy shows good 
convergence and approaches the GFMC result, in particular for the hQ = 22 MeV calcula- 
tion. The ground-state point-nucleon rms radius is almost HQ independent and converged. 
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It agrees with the GFMC value of 1.485 fm. 

We note that results on the first excited + state obtained using the resonating group 
method were reported in Ref. [171]. The Bonn potential employed in that work gives very 



similar results for the ground state as those obtained using the Argonne V8'. It is, therefore, 
reasonable to make a comparison for the excited state results. The first excited + state 
energy reported in Ref. |l7j was -6.42 MeV, which is about 10% below our result of iV max = 14 
and HQ = 16 MeV calculation. The reported rms radius, 3.02 fm, is slightly above our 
calculation. 

In order to further compare the convergence and the Q dependence of the results obtained 
with two- and three-body effective interactions, we present a similar calculation as that of 
Fig. ^|, obtained using the Argonne V8' NN potential and a larger Q difference, in Fig. |9|. The 
full lines correspond to the three-body effective interaction calculations, also shown in Fig. 
|7|, while the dashed lines connect the two-body effective interaction results. Two HO energies 
of HQ =16 and 22 MeV were used. The dotted line represents the GFMC result. Again we 
observe a better stability of the three-body effective-interaction results with respect to the 
model-space size changes, a smaller Q dependence, and a faster convergence, in paticular for 
the ground-state. In Table | we show the absolute value of the ground-state energy differences 
obtained in the calculations with HO energies of hQ = 16 MeV and HQ = 22 MeV for both 
the two-body and the three-body effective-interaction calculations in different model spaces. 
We can see that the differences obtained with the three-body effective interaction are almost 
two times smaller in model spaces with iV max = 6 — 10. The differences decrease with the 
enlargement of the model space for iV max > 8. We note that by construction the present two- 
and three-body effective-interaction calculations would become identical for iV max = 28. 

In Table || we present a summary of our results obtained in the largest model spaces used 
in the present study, e.g., iV max = 14 for the positive-parity states and iV max = 13 for the 
negative-parity states. The positive-parity state results were obtained using the three-body 
effective interaction. For the Argonne V8' NN potential calculations we also include the 
GFMC ground-state results for a comparison. We note that the Faddeev-Yakubovski 
equation solution gives -25.03 MeV || for the Nijmegen NN potential H35H , which gives com- 
parable results to the Reid93 NN potential for the three-nucleon problem. The experimental 
binding energy of 4 He is 28.296 MeV. The discrepancy between the experimental and cal- 
culated values are usually attributed to the real three-nucleon forces that were not taken 
into account either in our calculation or in the other calculations, which we discussed. We 
note that the difference in the binding energies obtained using the V8' and the Reid93 NN 
potentials is mainly due to the Coulomb interaction included in an isospin-invariant manner 
only in the calculations with the Reid93 NN potential. 



B. Charge form factors 

A sensitive test of the wave-functions obtained in our calculations is the evaluation of 
charge form factors. Using the formalism of Ref. [[I(J, we calculated the charge EM and 
strangeness form factors in the impulse approximation. The one-body contribution to the 
charge operator is given by Eq. (15) in Ref. e.g., 



13 



1 ^f 4°V) ., 



G ( e\t) 



2t^^-M, (25) 



where r 



2-j 



Li. is the fc-th nucleon orbital momentum, G$(t 



and I 



r) are the 
refers to 



one-body electric and magnetic form factors, respectively. The superscript (a 
(T = 0) for isoscalar EM form factor or to (s) for the strangeness form factor. We use the 



parametrization of the one-body form factors as discussed in Ref. |36j 



G%\r 



T 



(r 



G 



D, 



fJ>pG v (r 



VnTG°{T)t n {T) 



Di 



fi n G 
Ps tG°(t)Ut) 
» s G°(t) , 



(26a) 
(26b) 
(26c) 
(26d) 
(26e) 
(26f) 



with 



GviV 

tin 



:i + a£t) 

1 + A n r) 



(1 + A 



r 



(27a) 
(27b) 
(27c) 



The isoscalar EM form factor is given by G E 



\ i - \ [G^m + m\ i an d f° r the parameters 



appearing in Eqs. ( P6j , p7| ) , one has numerically fi p = 2.79, /x n = —1.91, Ay = 4.97, and 
A„ = 5.6. Following Ref. J36[, we also set the strangeness radius p s = —2.0 and A^ = A n . 
Limits on these parameters are to be determined in the experiments at the Thomas Jefferson 
Accelerator Facility (TJNAF). The first strangeness magnetic-moment measurement was 
reported recently [0 and an experimental value fi s = +0.23, obtained with a large error. 
We use this value in our calculations. 

Our charge form factor calculations are presented in Figs. 



10 -O. The charge form 



factors given in the figures were calculated using the one-body operator 



as 



2 v /?(/,0+0|il4 



r oo (?) ^ \i, + 0) . We show only results obtained with the Argonne V8' NN 
potential; the Reid93 NN potential gives almost identical results for the charge form factors, 
when the same HO energy HQ is employed. Our calculated elastic EM charge form factor 



is given in Fig. [10] together with the inelastic EM charge form factor corresponding to the 
transition to the first excited + state. These results were obtained using the HO energy 
hQ = 22 MeV and the three-body effective interaction in the iV max = 14 model space. In 
this calculation we obtained the best description of the ground state. The calculation of 
the elastic charge form factor in the impulse approximation can be directly compared to 



that presented in Fig. 2 of Ref. [36f], performed using Variational Monte Carlo (VMC) wave- 
functions and the Argonne V14 NN potential. There, the minimum was obtained at q ~ 3.55 
fm _1 , while the experimental minimum is at q w 3.2 fm -1 . The difference can be explained 
with the help of meson-exchange-current contributions. The elastic charge EM form factor 
obtained in our calculation compares well with that obtained by the VMC wave-functions. 
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It is shifted further to higher q, namely, we get the minimum at q ~ 3.75 fm _1 . We note that 
a second minimum appears in our calculated elastic charge form factor at q ~ 7.25 fm _1 . 
The second minimum at a similar position was found in the VMC calculations presented in 
Ref . |3£| . To examine the form factor dependence on hQ, we repeated these calculations for 



different choices of the HO energy. In Fig. [TTI we show the result obtained with Ml = 19 



MeV. All other characteristics are the same as in the calculation of Fig. [H]. The minimum 
here is shifted further to higher q, we have it at q ~ 3.85 fm -1 . The difference between 
the two results is rather small but still it shows that our calculation is not completely 
converged and, in particular, the description of the high transferred-momentum part of the 
form factors requires the use of even larger model spaces than we employed. We note that 
the inelastic form factor has a stronger Q dependence than the elastic form factor. As 
discussed in the previous subsection, convergence of the excited state has not been achieved 
in our calculations within the model spaces employed. Therefore, our calculated inelastic 
form factors must be taken with some degree of caution. Let us remark that, in addition to 
the transition form factor, we also computed the form factor of the first excited + state. 
That form factor was also evaluated in Ref. [|lj]] in the resonating group method approach 
using the Bonn potential. Similarly as in that work, our calculated 0^ form factor is almost 
an order of magnitude smaller than the ground-state form factor for a wide range of q. The 
first minimum is shifted in our calculation to larger q, more or less to the position of the 
transition form factor minimum, and the second minimum is shifted to smaller q compared 
to our ground-state form factor. 

Our calculated elastic strangeness form factor together with the inelastic EM charge form 
factor corresponding to the transition to the first excited + state are shown in Fig. [12. 
These results were obtained using the same wave-functions as those used for calculations 
presented in Fig. [It], namely we had hfl = 22 MeV, the model-space size characterized 
by N iaax = 14 and the three-body effective interaction was employed. The elastic form 
factor can be compared with the impulse approximation VMC result of Fig. 3 in Ref. |3T> 



Similarly, as for the EM elastic form factor, our calculation compares well with VMC result. 



We note, however, the different value of strangeness magnetic moment used in Ref. |36 
= -0.2). 



Finally, in Fig. 13 we present the ratio of the EM and strangeness form factors from 



Figs. [10] and [12]. The ratio of the elastic charge form factors is particularly interesting, as 
it can be experimentaly obtained from the measurement of the parity-violating left-right 
asymmetry for scattering of polarized electrons from a 4 He target. Experiments of this type 
are now under preparation at TJNAF. 



IV. CONCLUSIONS 

In the present study we used equations for Faddeev amplitudes, antisymmetrized for 
three nucleons, to solve the shell-model problem for the four-nucleon system. We performed 
calculations in larger model spaces, up to an HO excitation of 1AMI above the unperturbed 
ground-state, than in any other shell-model study so far. The main motivation for the 
present work was to test the shell-model approach and the effective interactions that we 
want to apply to more complex systems, e.g., p-shell nuclei, in particular. The effective 
interactions that we employed were derived from realistic NN potentials, i.e., the Reid93 
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and the Argonne V8'. In addition to the two-body effective interactions, we also computed 
the three-body effective interactions and demonstrated that their use significantly improves 
the convergence of the results. 

Our calculations depend on the model-space size and on the HO frequency Q. The 
effective interactions were constructed in such a way that in the large model-space limit the 
effective interactions approach the bare NN interaction. Thus our results should converge 
to the exact solutions. The dependence on the model-space size and Q was investigated. 
We found quite different behavior of the ground state and the first excited + state. Our 
ground-state energy and point-nucleon radius results begin to converge and are close to or 
in agreement with those obtained by the GFMC method. For the first excited + state, our 
results, the point-nucleon radius, in particular, show large model-space and fl dependence. 
This implies that significantly larger model spaces would still be needed in order to obtain 
the exact solutions. The nature of the state is discussed in the literature []i~4",16|. The 



Coulomb interaction plays an important role in the description of this state. In the present 
calculations we did not include the isospin breaking. Our formalism is quite general, however, 
and allows the use of interactions that break the isospin symmetry. On the other hand, the 
calculated properties of the 0^ state show better convergence behavior. In the model spaces 
studied, we obtained lower excitation energy of the 0J~ state than of the state, contrary 
to experiment. The extrapolation of the model-space dependence of these two energies to 
larger model spaces shows, however, that the correct ordering of the states will be obtained. 
Apparently, the state is associated with a radial excitation and, thus, it is more sensitive 
to the HO basis used in the expansion. 

A sensitive test of our calculated wave functions is the computation of the charge EM and 
strangeness form factors. Our impulse-approximation results show little dependence on the 
NN potential and our best results are close to the corresponding form factors obtained using 
the VMC wave functions and the Argonne V14 NN potential. In particular, we observe both 
the first and the second minima in the elastic charge form factor in positions close to those 
obtained using the VMC calculations. In addition to the elastic charge form factors, we also 
evaluated the form factors for the transition to the O2" state in the impulse approximation. 

In general, the energy scales of the bound 0s nucleons are significantly different from the 
scattering energies of the resonances. This difference can only be accounted for with a large 
N max in our approach. Consequently, the results for the excited states and the transition 
form factors obtained within the limited model spaces of the present work should be taken 
with some caution. In the future we would like to apply the formalism, discussed in the 
present paper, to a more extensive study of the negative-parity states of 4 He. In particular, 
it is desirable to use still larger model spaces to investigate the excited-state convergence 
properties. 

The most important result of the present work is, however, the successful use of the three- 
body effective interaction. This three-body effective interaction can be computed for more 
complex nuclei as well and, in principle, used, after a transformation to an appropriate three- 
nucleon basis, in standard shell-model calculations. A more practical approach, however, 
is to make use of the three-body effective interaction knowledge for the renormalization of 
the two-body effective interaction. Work in this direction is under way. In addition, the 
present formalism may be used to compute the four-body effective interaction for nuclei with 
A > 4. We plan to extend the shell-model Faddeev-like approach that we have successfully 
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applied to three- and four-nucleon systems to systems with more than four nucleons, using 
also a formalism of equations for components with lower degree of antisymmetry than the 
full wave-function developed in Ref. 
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FIGURES 



FIG. 1. The dependence of the ground-state and the first-excited + state energies, in MeV, 
on the maximal number of HO excitations allowed in the model space. The two-body effective 
interaction utilized was derived from the Reid93 NN potential. Results for h£l = 14, 17, 19 and 22 
MeV are presented. 

FIG. 2. The dependence of the point-nucleon rms radius of the ground state and the first-excited 
+ state, in fm, on the maximal number of HO excitations allowed in the model space. The 
two-body effective interaction utilized was derived from the Reid93 NN potential. Results for 
hfl = 14, 17, 19 and 22 MeV are presented. 

FIG. 3. The dependence of the ground-state and the first-excited + state energies, in MeV, 
on the maximal number of HO excitations allowed in the model space. Results obtained using the 
two-body (dashed line) and three-body (full line) effective interaction derived from the Reid93 NN 
potential are compared. Harmonic-oscillator energies of HQ = 17 and 19 MeV were used. 

FIG. 4. The dependence of the point-nucleon rms radius of the ground state and the first-excited 
+ state, in fm, on the maximal number of HO excitations allowed in the model space. Results 
obtained using the two-body (dashed line) and three-body (full line) effective interaction derived 
from the Reid93 NN potential are compared. Harmonic-oscillator energies of HQ = 17 and 19 MeV 
were used. 

FIG. 5. The dependence of the ground state, the first-excited + state and the first-excited 
CPO state energies, in MeV, on the maximal number of HO excitations allowed in the model space. 
For the positive-parity states the results were obtained using the three-body effective interaction 
derived from the Reid93 NN potential. Energies of the 0~0 state were calculated using a two-body 
effective interaction derived from the Reid93 NN potential. Harmonic-ocillator energies of = 17 
and 19 MeV were used. 

FIG. 6. The dependence of the point-nucleon rms radius of the ground-state, the first-excited 
+ state and the first-excited 0~0 state, in fm, on the maximal number of HO excitations allowed 
in the model space. For the positive-parity states the results were obtained using the three-body ef- 
fective interaction derived from the Reid93 NN potential. Energies of the CPO state were calculated 
using a two-body effective interaction derived from the Reid93 NN potential. Harmonic-oscillator 
energies of h£l = 17 and 19 MeV were used. 

FIG. 7. The dependence of the ground state and the first-excited + state energies, in MeV, 
on the maximal number of HO excitations allowed in the model space. The three-body effective 
interaction derived from the Argonne V8' NN potential was used. Results for fiQ = 16, 19 and 22 
MeV are presented. The dotted line represents the result -25.92 MeV of the GFMC calculation 

0- 
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FIG. 8. The dependence of the Point-nucleon rms radius of the ground state and the 
first-excited + state, in fm, on the maximal number of HO excitations allowed in the model 
space. The three-body effective interaction derived from the Argonne V8' NN potential was used. 
Results for hQ = 16, 19 and 22 MeV are presented. 

FIG. 9. The dependence of the ground-state and the first-excited + state energies, in MeV, 
on the maximal number of HO excitations allowed in the model space. Results obtained using the 
two-body (dashed line) and three-body (full line) effective interaction derived from the Argonne 
V8' NN potential are compared. Harmonic-oscillator energies of HQ = 16 and 22 MeV were used. 
The dotted line represents the result -25.92 MeV of the GFMC calculation Q. 

FIG. 10. The elastic EM charge form factor (full line) and the EM charge form factor corre- 
sponding to the transition to the first-excited + state (dotted line) calculated in the impulse 
approximation using the three-body effective interaction derived from the Argonne V8' NN poten- 
tial in the N mBX = 14 model space and HQ = 22 MeV. 

FIG. 11. The elastic EM charge form factor (full line) and the EM charge form factor corre- 
sponding to the transition to the first-excited + state (dotted line) calculated in the impulse 
approximation using the three-body effective interaction derived from the Argonne V8' NN poten- 
tial in the -/V max = 14 model space and hfl = 19 MeV. 

FIG. 12. The elastic strangeness charge form factor (full line) and the strangeness charge form 
factor corresponding to the transition to the first-excited + state (dotted line) calculated in the 
impulse approximation using the three-body effective interaction derived from the Argonne V8' 
NN potential in the N max = 14 model space and fi£l = 22 MeV. Values of the strangeness radius 
p s = —2.0 and the strangeness magnetic moment (j, s = 0.23 were employed. 

FIG. 13. The ratio of the elastic EM and strangeness charge form factors (full line) and the EM 
and strangeness charge form factor corresponding to the transition to the first-excited + state 
(dotted line) calculated in the impulse approximation using the three-body effective interaction 
derived from the Argonne V8' NN potential in the N max = 14 model space and HQ = 22 MeV. 
Values of the strangeness radius p s = —2.0 and the strangeness magnetic moment p s = 0.23 were 
employed. 
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TABLES 



■'"max 


6 


8 


10 


12 


14 




1.311 


1.466 


1.265 


1.037 


0.834 


\AE 3cS \ 


0.778 


0.782 


0.676 


0.601 


0.550 



TABLE I. Absolute value of the ground-state energy differences obtained in the calculations 
with HO energies of HQ = 16 MeV and HQ = 22 MeV with the two-body (second row) and the 
three-body (third row) effective interactions in different model spaces. The effective interactions 
were derived from the Argonne V8' NN potential. The corresponding energy dependence is shown 
in Fig. g. 



Argonne V8' NN potential 



State 


Variable 


m 


= 16 MeV 


nn 


= 19 MeV 


fin = 22 MeV GFMC 


Of 


E [MeV] 




-26.62 




-26.30 


-26.07 -25.92(8) 




vV) [fm] 




1.481 




1.485 


1.485 1.485(10) 


o 2 + 


E [MeV] 




-5.77 




-4.89 


-3.93 




E x [MeV] 




20.86 




21.42 


22.14 




J(r' 2 } [fm] 




2.906 




2.777 


2.658 


or 


E [MeV] 




-6.70 




-6.17 


-5.59 




E x [MeV] 




19.93 




20.14 


20.48 




V(r 2 ) [fm] 




2.349 




2.263 


2.186 


Reid93 NN potential 


State 


Variable 


nn 


= 17 MeV 


m 


= 19 MeV 




of 


E [MeV] 




-25.69 




-25.47 






vV) M 




1.487 




1.489 




0+ 


E [MeV] 




-5.00 




-4.39 






E x [MeV] 




20.69 




21.08 






vV) [fm] 




2.873 




2.787 




or 


E [MeV] 




-5.91 




-5.54 






E x [MeV] 




19.78 




19.93 






V(r' 2 ) [H 




2.339 




2.281 





TABLE II. Results for the ground-state and excited state energies and point-nucleon rms radii, 
as well as the excitation energies (E x ), obtained in the largest model spaces used in the present 
study, -/V max = 14, (13) for the positive- (negative-)parity states, respectively, are presented. All 
the states have isospin T = 0. The positive-parity-state calculations were performed using the 
three-body effective interaction. Results for different HO energies are given in separate columns. 



The GFMC ground-state results [34| are shown for comparison. 
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-16 
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-28 




-32 



19 +2 



17 x 

0"i 



Eeid93 



19 a + l 



3 4 5 6 7 8 9 10 11 12 13 14 15 16 

Nmax 



3.00 



^ 2.50 



5, 

§ 2.00 
* 1.50 



1.00 
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Nmax 



